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The nonlinear frequency shift coefficient, which represents the strength of the transformation of 
amplitude fluctuations into phase fluctuations of an oscillator, is measured for MgO-based spin- 
torque oscillators by analyzing the current dependence of the power spectrum. We have observed 
that linewidth against inverse normalized power plots show linear behavior below and above the 
oscillation threshold as predicted by the analytical theories for spin-torque oscillators. The mag- 
nitude of the coefficient is determined from the ratio of the linear slopes. Small magnitude of the 
coefficient (~ 3) has been obtained for the device exhibiting narrow linewidth (~ 10 MHz) at high 
bias current. 



Spin-torque oscillators (STOs) emit a microwave sig- 
nal. The signal originates in magnetization oscillations 
excited by bias dc current in a magnetoresistive (MR) 
device^. In recent years, extensive studies have been 
carried out on STOs because they are a promising can- 
didate for an on-chip microwave oscillator^!. One of 
the important properties of an STO is its frequency 
nonlincarity, i.e., a frequency depends on an oscillation 
amplitude 5 -^. Due to the property, the frequency of STO 
is tunable only by changing bias dc current, which is 
generally considered to be an advantage for applications. 
The nonlincarity is, however, a disadvantage of STO 
when thermal fluctuations are taken into account. Ac- 
cording to the analytical theory of Kim et al. 7 , amplitude 
fluctuations are transformed into phase fluctuations be- 
cause of the nonlinearity, resulting in spectrum linewidth 
broadening. The linewidth is a measure of the phase sta- 
bility of oscillation and it is preferable that it be narrow. 
Estimating quantitatively the nonlinearity, which deter- 
mines the device performances, is therefore a key sub- 
ject for further developments of STOs. According to the 
recent theories^&^ i 10 ' 11 ' 12 , the quantity representing the 
degree of the nonlincarity is the normalized dimension- 
less nonlinear frequency shift coefficient v (regarding the 
definition of the coefficient, see, e.g., Eq. (7) of Ref . Ho or 
Eq. (9) of Ref. [Tlh . In this letter, we report experimental 
estimations of the coefficient v in MgO-bascd STOs near 
threshold, which have not been directly addressed by pre- 
vious experiments. Following some theoretical remarks, 
experimental results are shown. 

It is theoretically known that the coefficient v has 
various values depending on magnetic environments and 
dampin g 10 ' 11 ! 12 . In general, the coefficient v also depends 
on bias curren t 10 ' 12 . Considering that the variation of v 
with bias current is small in the range near threshold, we 
treat v as a constant independent of bias current. 

First, we demonstrate numerically how large the coeffi- 
cient v is in typical STOs. A calculation example of v at 
the threshold in a planar device with uniaxial anisotropy 
is shown in Fig. [TJ representing the dependence on the 
magnitude and the angle of an in-plane external mag- 
netic field. The calculation is performed by the method 
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FIG. 1: (Color online) Dependence of v in the planar device 
with uniaxial anisotropy (i-direction) on the magnitude H 
and the in-plane angle 4>h of in-plane external magnetic field. 
The parameters chosen are an uniaxial anisotropy field Hk = 
200 Oe, a demagnetizing field 4ttM s = 8000 Oe, the Gilbert 
damping qg = 0.01, and the nonlinearity of damping^ qi — 
3. In the notation used here, positive and negative v denotes 
the red and blue frequency shift, respectively. 



described in Ref. [ll| which is based on the macrospin 
model. We have used the typical STO parameters shown 
in the figure caption. In the wide external field region, \v\ 
is much larger than unity {\v\ ~ 5-50). On the black line 
shown in Fig.[TJ the frequency nonlincarity vanishes (the 
nonlincarity due to the demagnetizing effect cancels out 
that due to the in-plane anisotropy), where remarkable 
reduction of linewidth is expecte d 12 ' 14 ' 15 . 

To estimate the value of v from experiment, we 
have used the spectrum analysis method based on the 
theories of spectrum linewidth of STO under thermal 
flue t uat ionsZ&^ i 1 ' 1 1 ' 1 2 . The method is similar to the one 
often used to measure the linewidth enhancement fac- 
tor (a-factor) in lasers 1 ^. According to the theories, the 
linewidth shows asymptotic behavior below and above 
threshold. In the region below threshold (the thermal 
activated oscillation region), the linewidth is given by 



A/< = A/ x (k B T/E osci ). 



(1) 



In the region above threshold (the current-induced oscil- 
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FIG. 2: (Color online) Linewidth and inverse power near 
threshold for several values of \v\. The units of linewidth 
and power are appropriately scaled. 

lation region) , the linewidth is given by 

Af> = Af Q x (k B T/E osci ) x (1 + is 2 )/2, (2) 

in which the additional phase diffusion due to amplitude 
fluctuations is expressed by the factor of v 2 . In Eqs. ([I]) 
and ©, A/o is the damping at thermal equilibrium ex- 
pressed in linewidth, k^T is the thermal energy and -Eosci 
is the average magnetization oscillating energy of a cer- 
tain oscillation mode. The energy E osc i becomes large 
as the bias current / increases. Since the voltage oscilla- 
tion signal from STO originates in the MR effect, -E^sci 
is proportional to a 'normalized' power P/I 2 , where P is 
the power of voltage oscillation signal from STO. There- 
fore, Eqs. ((T|) and © both indicate that the linewidth 
is proportional to the inverse normalized power in the 
two asymptotic regions, i.e., A/ < - > oc (P/I 2 )' 1 . Ac- 
cordingly, drawing a linewidth versus inverse normalized 
power plot, we can extract v from the following relation; 

3> / S< = {l + v 2 )/2, (3) 

where s< and s> are the slopes of asymptotes below 
and above threshold, respectively. To visualize Eq. ((3]), 
we show the linewidth and inverse power calculated nu- 
merically in Fig. [21 The calculation is based on the 
Fokkcr-Planck equation (FPE) corresponding to a noisy 
nonlinear auto-oscillator^ii. We have used the single- 
Lorcntzian approximation for a power spectrum. In the 
approximation, the linewidth is given by the real part of 
the first eigenvalue in the eigenmode expansion method 
for FPE (ReA M in the notation of Ref. H). In Fig. H 
linewidths show asymptotic behavior (dotted lines) rep- 
resented by Eqs. (Q]) and ^ in the region below and 
above threshold, and the relation between the slopes of 
asymptotes (Eq. Q) holds. 

Measurements are performed on MgO-bascd planar 
magnetic tunnel junction devices at room tempera- 
ture. Tunnel junctions are composed of IrMn(10)/ 
PL/MgO(1.05)/FL, in which the free layer (FL) is a 
CoFeB(3) layer and the pinned layer (PL) is a CoFe(4)/ 
Ru(0.95)/CoFcB(4) synthetic antifcrromagnct trilaycr. 
The junctions have been fabricated by sputter deposi- 
tion with annealing for one hour in a high magnetic field 
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FIG. 3: (Color online) PSD of the voltage oscillation signal 
in the setup (I) for bias currents from 0.4 to 0.7 mA. 



and at the temperature of 300° C. Nanopillar devices 
are patterned using electron-beam lithography and ion 
milling. On the same wafer, we find samples with vari- 
ous TMR ratio in the range 10-90%. We have observed 
that lower TMR ratio samples tend to show narrower 
spectral linewidth oscillations. This tendency in MgO- 
based samples is similar to that reported in Ref. [H. We 
consider that the difference of TMR ratio among sam- 
ples may result from localized defects in the MgO barrier 
as speculated by the authors of Ref. In this letter, 
we show the results of one typical low TMR ratio sam- 
ple with a 200 x 120 ran 2 elliptical shape. From R-H 
(resistance-field) characteristic measurements of the sam- 
ple for several field angles, we find that the resistance is 
235 (5.6 n • fim 2 ) for parallel configuration, TMR ratio 
is 15.6 %, and an in-plane anisotropy field Hk is about 
200 Oc. The pinned layer magnetization is along the long 
axis of the ellipse without an applied field. 

Oscillation properties are measured by a spectrum an- 
alyzer. We apply dc current / to the sample with an- 
tiparallcl magnetization configuration. Positive / corre- 
sponds to electrons flowing from the pinned layer to the 
free layer. Oscillation properties vary sensitively with 
the field applied in the sample plane: the magnitude H 
of 200-700 Oe and the angle (j) H of 0-90°. When the angle 
<t> H oi 40° ± 5° and the magnitude H ~ 500 ± 50 Oe, dis- 
tinct oscillation peaks are observed. In the other applied 
fields, distinct peaks are not observed at any bias current 
(/ < 0.7 mA). Below, we show the results for the two 
setups of the applied fields in which narrow linewidths of 
- 10 MHz are observed: (I) H = 490 Oc and <j> H = 42.4°, 
and (II) H = 516 Oc and cj) H = 39.7°. 

Power spectral density (PSD) of the voltage oscillation 
signal in the setup (I) is shown in Fig. [31 The character- 
istic oscillation mode is observed around 7.2 GHz for the 
bias current / > 0.42 mA. This mode grows steeply as 
the bias current increases and exhibits narrow linewidth 
in the high-current region. The PSD in the setup (II) 
is similar to that in the setup (I). The dependence of 
the frequency, the full width at half maximum (FWHM) 
linewidth, and the normalized power of the mode on the 
bias current for the setups (I) and (II) is shown in Fig. 2) 
The data are obtained by Lorentzian fits to the spectral 
peaks of the characteristic oscillation modes. The bro- 
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FIG. 4: Frequency ((a) and (d)), FWHM linewidth A/ ((b) 
and (e)), and normalized power P/I 2 ((c) and (f)) vs. bias 
current I for the setups (I) and (II). In (b), A/ ~ 14 MHz at 
0.58 mA. In (e), A/ ~ 8.6 MHz at 0.56 mA. 
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FIG. 5: (a) and (b) A/ vs. {P/I 2 )' 1 -plots for the setups (I) 
and (II). 



ken lines in Fig. |4] denote the threshold currents ith- The 
normalized power below threshold depends on the bias 
current in the way that (P/I 2 )^ 1 oc (7 t h — I)& By using 
the expression, we find that ith — 0.46 mA for the setup 
(I) and 7 th - 0.44 mA for the setup (II). In Figs. H (a) 



and (d), the frequency shift with the bias current is much 
smaller than that observed in previous experiments for 
MgO-bascd STO o 18 i 19 , and we cannot judge whether the 
shift is red or blue. In Figs. 2] (c) and (f), the normalized 
power grows steeply up to the current I ~ 0.52 mA from 
the threshold current. For I > 0.52 mA, the normalized 
power is being saturated. 

From the data for linewidth A/ and normalized power 
P/I 2 shown in Figs.|H we have obtained the A/ vs. (P/ 
i 2 ) _1 -plots shown in Figs. [3 We have used the data 
for the current I < 0.52 mA. As predicted by the the- 
ories, we observe that linewidth against inverse normal- 
ized power plots show linear behavior below and above 
the threshold. The appearance of the linear behavior 
above the threshold supports the consideration that the 
variation of v with bias current is small in the range near 
threshold. By estimating the values of linear slopes s<,> 
and using Eq. ([3]), we find the coefficient \v\ ~ 2.9 for the 
setup (I) and \v\ ~ 2.6 for the setup (II). Considering 
the range of values expected by the calculation as shown 
in Fig. [H these values of the coefficient (\v\ ~ 3) are 
comparatively small. The smallncss is consistent with 
the flatness of frequency (Figs. [4] (a) and (d)) and the 
exhibition of narrow linewidths of ~ 10 MHz in the high- 
current region (Figs.Q](b) and (e)). According to Fig.[TJ 
the coefficient \v\ is expected to have small values when 
an in-plane field is applied to the hard axis, i.e., 4>h ~ 80- 
90°. Our results measured in the angle of 4>h ~ 40° are 
inconsistent with the expectation. We consider that the 
mode with 7.2 GHz peak may be a non-uniform mode 
and so the effective directions and structures of in-plane 
anisotropy for the magnetization composing the mode are 
complicated. Micromagnetic study to clarify the oscilla- 
tion mode is now in progress. 

In summary, we estimated the coefficient v in MgO- 
based STOs which is a measure of the transformation 
of amplitude fluctuations into phase fluctuations of an 
oscillator. For the device exhibiting narrow linewidth 
(A/ ~ 10 MHz) in a high-current region, small value of 
the coefficient [\v\ ~ 3) was obtained. 
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